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Abstract 

In this paper the dynamics of Pierce diode with overcritical current 
under the influence of delay feedback is investigated. The system 
without feedback demonstrates complex behaviour including chaotic 
regimes. The possibility of oscillation regime control depending on the 
delay feedback parameter values is shown. Also the paper describes 
construction of a finite-dimensional model of electron beam behaviour, 
which is based on the Galerkin approximation by linear modes expan- 
sion. The dynamics of the model is close to the one given by the 
distributed model. 

1 Introduction 

Pierce diode is one of the classical models of plasma microwave electronics 
[1-9]. This distributed model though is rather simple, demonstrates many 
features of the electron beam dynamics in different real electron devices. 
The model consists of two infinite parallel plains pierced by monoenergetic 
electron beam. The region between two plains is evenly filled by neutralizing 
stationary ions, which density is equal to the non-perturbed beam electron 
density. The only controlling parameter is named Pierce parameter 

a = ujpL/v 0j 
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where oo p is the electron beam plasma frequency, L is the distance between 
diode plains, Vq is the non-perturbed electron velocity. It was already shown 
that with a > 7i so-named Pierce instability develops in the system and the 
virtual cathode is formed in the electron beam [2,3]. At the same time in a 
narrow range of Pierce parameter near 3tt the increase of the instability is 
suppressed by the non-linearity and in the electron beam the regime with- 
out reflection realizes. In this case the system may be described by fluid 
equations. It was also shown that in a narrow range of Pierce parameter the 
system may represent chaotic dynamics [4-9]. 

Recently the possibilities of oscillation control in finite-dimensional dy- 
namical systems are explored in detail [10-14]. At the same time the problem 
of distributed systems regime control causes great interest. For example it 
may be realized by adding to a system a controlled delay feedback [15,20]. 
The influence of the delay feedback on the dynamics of the beam with over- 
critical current in the regime of virtual cathode oscillation have been already 
investigated and the influence of the feedback on the characteristics of gen- 
eration have been shown [16]. The problem of chaos control in Pierce diode 
has also attracted the attention. In particular, the early work [17, 18] con- 
sider the possibility of stabilizating chaotic dynamics in Pierce diode in the 
regime without electrons' reflection and in the regime with virtual cathode 
with the help of E. Ott, C. Grebogy, J. Yorke [10] method. In our work [19] 
we analyse the possibility of stabilizating of chaotic dynamics of the fluid 
model of Pierce diode using the continuous delay feedback [11]. 

But the practical realization of such schemes in microwave devices come 
across several difficulties. So the investigation of chaotic dynamics of dis- 
tributed active media with external delay feedback causes great interest. 

So it seems interesting to examine the delay feedback influence on the 
dynamics of the system without virtual cathode - the fluid model of Pierce 
diode - because it represents all the classical regimes of a real distributed 
self-oscillation model. In this work the influence of the feedback parameters 
to the system characteristics is examined. The distributed chaotic system is 
analysed with the help of the numerical modelling of the original system of in 



2 



partial derivative non-linear equations and with the use of finite-dimensional 
model. 

The structure of the work is the follows. In section 2 we discuss the fluid 
model of the Pierce diode without feedback. In section 3 explored in detail 
the dynamics of this model with the added delay feedback depending on the 
value of delay time and feedback amplitude. The finite dimensional model 
of the investigated system using the Galerkin method is constructed and its 
behavior under the influence of the delay feedback is investigated in section 
4. The dynamics of finite-dimensional model is compared with the behaviour 
of the electron beam in the distributed system. 



2 Explored model 

The dynamics of Pierce diode processes in fluid electronics approximation is 
described by movement, continuity and Poisson equations: 

dv dv dip 
dt dx dx ' 

dp dp dv 

^ +v ix +p irx = > (2) 



d 2 lf 



a\p-l\ (3) 



dx 2 

where <p is the space charge field potential, p is the electron density, v is the 
electron beam velocity. 

The boundary conditions are: 

v(0,t)=v o , p(0,t)=p , <p(0,t)=<p(l,t) = 0. (4) 

The initial conditions are taken as a small perturbation of the space charge 
density near the homogeneous equilibrium state 

v(x) = v , p(x) = po, vi x ) = 

as p(x, 0) = psin27ra; where p 1. The equilibrium state becomes unstable 
then a > ii. In the equations the normalized values are used. 
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The delay feedback is brought in by modelling the potential difference 
between entrance and exit grids by the signal taken off from the interac- 
tion space in the point x = x^. As a control signal the oscillations of the 
space charge density p(xdf,t) is used. It can be interpreted as connecting a 
waveguide with a delay line to the interaction space, which is excited by the 
electron beam oscillations. Adding the delay feedback into the model leads 
to the changes in the right boundary conditions 

¥>(1, t) = f d f(t) = A(p(x df , t-d)-po). 

Here A is the delay feedback coefficient, characterizing the part of the oscil- 
lation power branched to the feedback delay line, d is the delay time value. 
Assuming that the development of the processes in our system begins at t — 
and when t < the space charge density is non-perturbed p(x,t) = po, the 
initial distribution of the delay feedback function is written as 

fdf(t)\te[-d,o] = 0. 

We have found out that the point of connection x^f does not influence 
upon the dynamics of the system. In this work value xjf is fixed as x,if = 0.2. 

Numerical solution of the equations (1) and (2) is found using explicit 
scheme with differences against flow. Poisson equation (3) is integrated using 
error vector propagation (EVP) method [21]. 

3 Delay feedback influence on the nonlinear 
dynamics of electron beam 

In Pierce diode without delay feedback when a decreases from 2.887T to 2.867T 
the behaviour of the electron beam changes from regular via period doubling 
cascade to weakly chaotic with neatly expressed time scale. With the further 
decrease of Pierce parameter the chaotic oscillations of the beam complicate 
essentially, the time scale disappears and spectral distribution complicates. 
We call this two types of chaotic behaviour ribbon and spiral chaos. All 
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the results described in this paper have been derived for a = 2.867T, i.e. a 
system without delay feedback must represent the "spiral chaos" oscillations. 
As quantitative characteristics of oscillation regime correlation dimension 
D [22] and highest Lyapunov index A [23] for the restored attractor are 
taken. This values do not change in the different points of interaction space 
D = 2.18 -J- 0.01, A = 0.16 -J- 0.04. In figure 1 the different regimes are 
represented on the A-d parameter plane with a = 2.867T. 
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Figure 1: Oscillation regimes reproduced on parameter plane (A, d) (a = 
2.86tt, x df = 0.2) 

For comparison the non-dimensional characteristic time of oscillations in 
the electron beam r = 4.06. The areas of n-periodical oscillations on the 
parameter plane are marked as T n When i <C 1 the system demonstrates 
chaotic oscillations identical to those without feedback. 

In figure 2 the phase-plane portraits, spectrums and time series for the 
"spiral" chaos regime (a), weakly developed chaos (b) and regular (period 1) 
oscillations (c) are represented. The phase-plane portraits are reconstructed 
by Takens delay method [24] from the time series characterizing the oscilla- 
tion of the space charge density in the fixed point of the distributed system. 
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Figure 2: Phase portraits, spectrums and time series of electron beam oscil- 
lations for the cases of "spiral" chaos regime (a), weakly developed chaos (b) 
and regular (period 1) (c) regimes 

Though system behaviour (a) is rather complex, in the frequency spectrum 
one can see the base frequency fo = 1/r = 0.25 and its second harmony 
2_/o • Analyzing the system attractor in phase-plane space one can see that 
near the instable state p = po there is a loop on which the motion of the 
phase point become slower. The other space of attractor is tightly filled by 
spiral phase pathes. When A increases the behaviour of the system may be 
different depending on the value of delay time d. If d > r/2 the complexity 
of electron beam dynamics increases with the increasing of A. The frequency 
spectrum and the phase-plane portrait become more complex too, the os- 
cillation amplitude enlarges. Further enlargement of A leads to essentially 
different dynamics of the system. The oscillation amplitude sharply increases 
and then in the electron beam the reflected electrons appear. In the system 
appears the virtual cathode. The electron beam behaviour is determined by 
two mechanisms - the Pierce instability and its limitation by the nonlinearity. 
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The delay feedback with parameter values d > 1.8 2 A > 0.02 destroys the 
limitation mechanism and leads to the increasing of instability and further 
to virtual cathode forming. In this case the fluid model becomes incorrect 
because equations (1), (2) describes the processes in the electron beam only 
without overtaking or reflection. The threshold value Aye depends on the 
value of d. In the case d < r/2 the complexity of the electron beam oscillation 
decreases with the increase of A. The noise base diminishes, with the further 
enlargement of A on the bifurcation map one can see periodic gaps. When 
A > 0.03 a transition from the chaotic to periodical dynamics via reverse 
doubling period cascade have place. The oscillation amplitude decreases and 
approach to non-perturbed value po- In a wide range of feedback parameters 
it is possible to suppress the chaotic dynamics of the electron beam and to 
establish the regular one. In figure 1 the areas of cycles 1-16 period and tran- 
sitions between different regimes depending on the changes of delay feedback 
parameters are shown. 

Now some words about the physical processes in the electron beam. As 
the explorations of the electron waves propagation show, dynamics of electron 
beam is mainly determined by the distance between the current state of 
the system and the homogeneous equilibrium state p(x) = p , v(x) = vo, 
(p(x) = <fo. The distance between the current state of the system and the 
homogeneous equilibrium state can be determined as 

S(t) = ^J (p(x, t) - po) 2 + (<p(x, t) - v? ) 2 + (v(x, t) - v f dx j . (5) 

Time-dependent changes of this value illustrates figure 3 for the cases of 
spiral chaos regime (a), weakly developed chaos (b) and regular (period 1) 
oscillations (c). 

In the first case the system in some time comes very close to the homo- 
geneous equilibrium state S ~ and the oscillation amplitude is very small. 
Then the mechanism of instability activates and the amplitude of the os- 
cillations increases until it's limited by the nonlinearity. Then the process 
repeats, but each time S and the space distributions of the values near the 
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Figure 3: The distance between the current state of the system and the 
homogeneous equilibrium state depending on time for the cases of "spiral" 
chaos regime (a), weakly developed chaos (b) and regular (period 1) (c) 
regimes 

equilibrium state are different, so the development of instability begins from 
another conditions and the dynamics of the system is irregular. The dynam- 
ics of the system can also be examined by considering non-linear energetic 
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functional [25]: 

1 1 

Wk = \j pv2 dx ~ h Wp = \j pvdx ^ 

o o 
This functions describe energy transitions between kinetic energy of the 
beam movement and potential energy of the space charge field. One can see 
that for the chaotic regimes the maximums of the functionals are larger than 
those for the regular processes, because of the larger degree of non-linearity 
in the chaotic regime. 
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Figure 4: Difference between kinetic energy and potential energy of the beam 
movement AW = Wk — W p for the cases of "spiral" chaos regime (solid line) 
and regular regime (dotted line) 

Figure 4 illustrates the time- dependence of the function AW = W k — W p 
characterizing the energy transition processes in the electron beam. In the 
periodic regime (the dotted line) the energy transitions are regular. The 
maximums of the functions shows the accumulation of the charge in the 
interaction space. In the chaotic regime (the solid line) in some moments the 
energy difference AW = 0. The system nears to the unstable equilibrium 
state. Then the wave movement energy increases again abruptly. Further 
the perturbation energy decreases and the system protractedly is situated 
near the equilibrium state. It can be physically explained as follows: in the 
electron beam the stationary wave of electron density appears. Its amplitude 
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increases abruptly near the exit grid. Then the electromagnetic field of this 
clot brakes the following electrons and the result is that much more electrons 
leave the interaction gap then enter it. The discussed values approach the 
equilibrium state and then the process repeats from the new starting point. 
The feedback destroys this mechanism. When the delay time d < r/2 the 
system cannot approach the equilibrium state because the feedback signal 
extenuates the stored charge in the exit region and promotes the acceleration 
of the beam in the instant time when AW is maximal. And vice versa when 
d < r/2 the feedback signal leads to the increasing of oscillation amplitude, 
the system dynamics become more complicate and with sufficiently large A 
the virtual cathode appears. 



4 Finite-dimensional model of electron beam 
dynamics 

In [2] the method of constructing a finite-dimensional model based on the 
Galerkin approximation by linear modes expansion is described. It has been 
shown by [26] that in the range of Pierce parameter variation a G (2tt, 3tt) in 
the system excites infinite number of modes which can be determined from 
the dispersion equation 

{expfjaw] [(zu 2 + 1) sin a + 2jzucosa\ + am 4 — 

-aw 2 - 2jvu} [w 2 - I)- 2 = 0, (7) 

where w = u/uo p . In the case of Pierce diode system the modes were deter- 
mined by Kuhn (1986). It have been shown that among the excited modes 
only three were damping rather slowly and containing the most part of the 
system energy. So for the description of the system dynamics it is enough to 
take into account only this three modes. For the different Pierce parameter 
values corresponding to different dynamical regimes the space distributions 
vary weakly, so we can suppose they are independent of a. The initial basis 
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for the finite-dimensional approximation is taken as 

3 

i=i 

3 

p = J2M*Xt), (8) 

i=i 

3 

i=i 

where Ri, Vi, $j are the space distributions modes for the Aj, on are the modes 
amplitudes. Substituting the trial solution (8) into system (l)-(3), written 
for weakly perturbed values, we derive the nullity vector \P = (\I/i, \I/2, ^3), 
which components can be written as 

tfi = difli + d 2 R 2 + a 3 R 3 + a x {R lx + V lx )+ 

+a 2 {R 2x + V 2x ) + a 3 {R 3x + V 3x ) + 
+a 2 1 (R 1 V 1 ) x + a 2 2 (R 2 V 2 ) x + a 2 3 (R 3 V 3 ) x + 
+a 1 a 2 (R 1 V 2 + R 2 V ± ) X + a 1 a 3 ( J R 1 V r 3 + i? 3 ^i)x+ 

+ a 2 a 3 (i? 2 \/ 3 + #3^2)*, (9) 

^ 2 = diVi + d 2 V2 + dgVg + + + 

+a 2 ($2. + V 2:c ) + a 3 ($ 3:c + V 3x )+ 
+0^1^ + a^V^ + a 2 3 V 3 V 3x + 
+a 1 a 2 (V 1 V 2x + V 2 V lx ) + a^V^ + V 3 V lx )+ 

+ a 2 a 3 (V 2 V 3x + V 3 V 2x ), (10) 
^ 3 = ai($i^ + a 2 R{) + a 2 ($ 2xx + a 2 R 2 ) + 

+ a 3 {<5> 3xx + a 2 R 3 ), (11) 

where (•)* = J^(-). 
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The internal product of functions is defined as 



(fxg)= fg dx. 



(12) 



Using Galerkin method we can find the unknown coefficients from the 
matrix equation 

/ / R ± V 1 $j. \ / *i \ \ 



V 



R 2 V 2 $2 
R 3 V 3 $ 3 



X 



0. 



(13) 



^2 

Carrying out elementary transformation and taking into account the equa- 
tions (9)— (11), we derive the matrix equation for the coefficients af. 



MA + BA + D = 0, 



(14) 



where vector A is composed from coefficients a^. The elements of matrixes 
M and B are derived as 

mi d = (Rj x Ri) + (Vj x Ri), 

bij = ((R ix + V ix ) x Rj) + ((V ix + $«) x Vj) + 
+(($ ixx + a 2 Ri) x $,). 
Matrix element D is derived from formula 

di = Yl a l i( R kVk) x Ri + V k V kx Vi\ + 
k 

+ mi akai ^ RkVi + R i v k)xRi+ 
k i, \+k 

+{y k v lx + vy^M . 

Resolving the equations (14) relatively dj, we derive the explicit equations: 
d\{t) = Z n ai + /i 2 a 2 + /i3a 3 + li±a\ + h<ja\+ 

+ ^16^3 + li7a 2 a 3 + hsaiai + Zi 9 aia 3 , (15) 
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a 2 {t) = Z 2 iOi + /22O2 + ^23^3 + ^24«l + h5<4+ 
+ /26°3 + ^27 a 2°3 + ^28°l a l + ^29°l a 3; 

03(0 = hlUl + ^320 2 + ^33^3 + kAa\ + ^35« 2 + 
+ ^36 a 3 + ^37 a 203 + ^38°l a l + ^39°l a 3- 



(17) 



(16) 



The coefficients Uj are derived from numerical solution of the implicit equa- 
tions and are independent of Pierce parameter variations. The non-linearities 
of the system are quadratic and appears because of the kinematic non- 
linearities and those presenting in continuity equation. The 1 and 2 modes 
are excited by the instability negative dissipation and its energy is inherited 
into the 3rd mode. As the numerical analysis shows the finite-dimensional 
demonstrates the same types of behaviour as the distributed one. With the 
decrease of a the system dynamics becomes more complex and further tran- 
sition to chaos via period doubling cascade takes place. The comparison 
of bifurcation diagrams for the distributed model and the finite-dimensional 
model is made in figure 5. 

Also in the system two variants of the chaotic regime similar to those 
in the distributed model are observed - the "spiral" chaos and the "band" 
chaos. The delay feedback is brought in by adding into the right part of the 
equations (7)- (9) the signal F d f(t — d) which is formed as: 



In figure 6 the bifurcation map in the parameters A - d for the finite- 
dimensional model is presented. The Pierce parameter a = 2.7747T. With 
this parameter value in the system without feedback the "spiral" chaotic 
oscillations are observed. In the map we can see that with the increase of 
feedback signal amplitude and d < r/2 (where r m 3.15 = l//o, fo - the base 
spectrum frequency) the chaotic dynamics is suppressed and regular regime 
is installed. 

When d > r/2 the oscillation amplitude sharply increases which is equal 
to system transition to virtual cathode forming regime, where the finite- 
dimensional equations become incorrect. Comparing figure 1 and figure 4 



F df (t) = A[ ai (t) + a 2 (t) + a 3 (t)} 



(18) 
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Figure 5: Bifurcation diagrams for the distributed model (a) and the finite- 
dimensional model (b) 
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Figure 6: A-d parameter plane for the finite-dimensional model with delay 
feedback, a = 2.86n 
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one can see that the finite-dimensional model gives a very good description 
of the processes taking place in the distributed system. 

5 Conclusion 

In our work the delay feedback influence on the electron beam dynamics in 
hydrodynamical and finite-dimensional models of Pierce diode is investigated. 
It is shown that with some feedback parameters' values the chaotic dynamics 
of the electron beam is suppressed and periodical regimes of different types 
may be installed. Physically it is connected with the changing of conditions 
of electron waves propagation. The practical interest for this phenomenon is 
caused by the ability of eliminating the undesirable parasitical and chaotical 
oscillations in some real systems where Pierce instability may appear (for 
example, in electron guns, beams of charged particles, etc.). 
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